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We show how dynamical heterogeneities in glass forming systems emerge as a consequence of 
the existence of dynamical constraints, and we offer an interpretation of the glass transition as an 
entropy crisis in trajectory space (space-time) rather than in configuration space. To illustrate 
our general ideas, we analyze the one dimensional (d = 1) Fredrickson-Andersen and East models. 
Dynamics of such dynamically constrained systems are shown to be isomorphic to the statics of d-\- 1 
dimensional dense mixtures of polydisperse non-interpenetrating domains. The domains coincide 
with arrested regions in trajectory space. 
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A glass forming system, like a supercooled liquid, ex- 
hibits a precipitous onset of slowness. As temperature 
is decreased in these systems, typically in a range of a 
few decades, relaxation times and viscosities increase by 
several orders of magnitude, eventually surpassing exper- 
imentally accessible times. For practical purposes, these 
systems effectively freeze at the glass transition temper- 
ature T g . For reviews see |J. Interestingly, this dynam- 
ical arrest carries no evident static structural signature 
of growing length scales. Rather, experiments and sim- 
ulations show that supercooled liquids are dynamically 
heterogeneous [[||| . Molecules in one region of the liquid 
translate or rotate several orders of magnitude faster or 
slower than those in a neighboring region. The spatial 
extent of these dynamical heterogeneities is mesoscopic, 
and the time scale of the slowest domains increases with 
decreasing temperature at least as fast as the relaxation 
time of the system. Such structural behavior seems be- 
yond description with homogeneous methods like mode 
coupling H and mean field theories ||, and it is widely 
neglected in analytical treatments (see, however, 
Nevertheless, we show here that for a broadly applicable 
mechanism of dynamical arrest, these heterogeneities are 
intrinsic to the nature of glass forming systems. 

Our central result is that dynamical heterogeneities are 
a manifestation of the existence of nontrivial structure in 
the trajectories of glassy systems. This structure associ- 
ated with dynamics is independent of any specific static 
properties. Instead, the nontrivial dynamical structure is 
a consequence of local dynamical rules that significantly 
restrict the size of accessible trajectory space. For exam- 
ple, consider a highly compressed (or supercooled) glass 
former. Atoms in most regions of space are jammed, 
making mobility possible in only a relatively low fraction 
of spatial regions. These rare regions are those that are 
already unjammed, or those that may be close in space to 
an unjammed region. In the evolution of such a system, 
one therefore expects a clustering of mobile regions and 
thus a mesoscopic demixing of mobile and static regions. 
Macroscopic demixing is not expected because dynamics 
should conserve a canonical distribution. This picture of 
the origin of dynamic heterogeneity is in accord with the 



idea that glassiness is not necessarily a consequence of ei- 
ther disorder or frustration in the static interactions but 
of the existence of effective constraints on the dynamics 
of the system 

The simplest microscopic models that illustrate this 
view are the Fredrickson-Andersen (FA) model ]l0| and 
the East model (TTJ . They consist of a chain of Ising 
spins <Ji = ±1 {i = l---iV), with trivial Hamiltonian 
H = oi, and single spin flip dynamics subject to local 
kinetic constraints. In the FA model, a spin can flip if 
either of its nearest neighbors is in the up state jHJ. In 
the East model, a spin can flip only if its nearest neigh- 
bor to the right is up |ll|] . The equilibrium behavior 
of both models is that of an uncorrelated spin system, 



with c = 1/ (1 



ol/T\ 



being the average concentration 



of up spins at temperature T. The competition between 
decreasing the energy and the need for facilitating spins 
leads to a glassy slow down at low temperatures. The 
relaxation times go as r ex e 3 / T in the FA model, and as 
the East model. These models therefore 
correspond to strong and a fragile glass formers, respec- 
tively. See [Q for a short review. 

Let us first establish that both the FA and East models 
display dynamical heterogeneities. Consider the coarse 
grained spins Si(t;At) = (Afp 1 Jo dt'oi(t + Fast 
and slow spins in the time window of width At at time 
t will correspond to low and high values of (t; At), re- 
spectively. The spatial distribution of sf(t; At) will de- 
termine the extent to which the system is dynamically 
heterogeneous. The heterogeneity made evident with 
this field depends on the coarse graining time. For At 
very short, only the trivial uncorrelated static structure 
is probed. The same is true for At much larger than the 
relaxation time where ergodicity is restored. For interme- 
diate values of At, however, we may expect to see spatial 
structure in the (t; At). I n Fig. 1 we show that this is 
indeed the case by plotting the structure factor, S(k), for 
the field sf(t; At) calculated from simulations of the FA 
model (left panel) and East model (right panel). We have 
used coarse graining times of about a fifth of the relax- 
ation time, and simulations were performed using contin- 
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uous time Monte Carlo ]l3| for systems of 10 5 spins av- 
eraged over 10 3 samples. For both models, the existence 
of a correlation length £a£ is clear. Moreover, grows, 
albeit slowly, with decreasing temperature and thus in- 
creasing relaxation time. For the FA model the structure 
factor decays approximately as S(k) ~ k~ 2 for large k. 
In the East model, at intermediate fc, the structure fac- 
tor goes as S(k) ~ ^- in 3/ in 2 ^ p or } ar g er momenta, the 
structure factor is oscillatory, a feature which becomes 
more pronounced with decreasing temperature. 

In order to understand these behaviors, we study the 
generating functional for the trajectories of these sys- 
tems. Consider a trajectory {<j t } = (<to, <ti, <?2, . . . , ar), 
where a t = (ait, ait, • • • , crjvt) denotes the configuration 
of the TV spin system at time t. Each trajectory has a 
probability P({at}), so that the total probability of going 
from <Jo to Br is the weighted sum of such trajectories, 
^ ,ar = E{a t } P ({^})- This partition sum, Z# 0i z r , 
can be written as |Ol 
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The first factor in the summand of (1) corresponds to 
the probability of a trajectory in the absence of dy- 
namical constraints. Its action is just that of N non- 
interacting ferromagnetic Ising chains in uniform fields, 

S o({a t }) = EiIlEt=0 1 ( J + VitVit+l + J- + -(Tit + 

J — + cTit+i + where J^ V£ = {ln[l - j5t (1 - c)] 

+ /i ln(l -jStc) + 1/ In [7 ft (1 - c)] + e ln(7 5t c)} /4, 
7 is a microscopic rate, and 5t the time step. The sec- 
ond factor in the summand of (1) embodies the kinetic 
constraint. Only trajectories which satisfy the condition 
hu = at all space-time points are allowed. If hu is 
chosen to be a positive number when the constraint is 
not satisfied, then Y\ it S(hit) = S(H) = e~ xn , where 
TL = Y^it hit and A — > 00. The infinite coupling con- 
stant means that the space of allowed trajectories is that 
of the ground states of TL. The respective forms of TL for 
the FA and East models are 

^FA = ^2 (1 _ VitVit+l) (1 — Gi+lt) (1 — ai-it) /8 , 
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^ (1 - a it a it+1 ) (1 - a i+lt ) /4 , 



so that the i-th term is zero unless the spin i flips and the 
facilitating spins are in the down state. The third factor 
in the summand of Eq. (1) accounts for preservation of 
norm and detailed balance. Trajectories which are al- 
lowed by the constraints have different probabilities than 
they would in the unconstrained case. In particular, 

AS FA = ( X - ^+1*) ( X - ^i-i*) (j-^t + J+) /4 , 

it 

A^East = ^ (1 ~ CTi+lt) (j-^it + J+^) /2 , 



where J± = J ±00 + ln[l - jSt (1 - c) e~ A ]/4 ± ln(l - 
7#£ce -A )/4. These expressions describe a proper dy- 
namics for any A. When A = 0, the expressions reduce 
to the trivial model of independent spins with uncon- 
strained dynamics. When A — > 00, they correspond to 
the FA and East models. 

TL and AS introduce competing spatial interactions in 
the space of trajectories. Those of TL are strong, and they 
are ferromagnetic in the sense that they favor the clus- 
tering of like spins. In contrast, those of A£ are weak 
and antiferromagnetic. Moreover, the scaling with dis- 
tance of interactions in TL is different than that for those 
in AS, so that nontrivial structure in the space of tra- 
jectories can be expected. In Fig. 2, we show samples of 
equilibrium trajectories for the unconstrained case (top), 
the FA (bottom left) and East models (bottom right), 
at T — 1.0. The difference between the constrained and 
unconstrained dynamics is striking. The trajectories in 
both the FA and East models display an extensive num- 
ber of domains of down spins. These domains are the 
origin of the dynamical heterogeneities. Spins within 
these domains do not change, so when coarse grained 
in time, they correspond to slow regions. To the extent 
that lowering temperature and thus energy in an atomic 
system coincides with decreasing facilitating spin concen- 
tration c, the structures observed in Fig. 2 coincide with 
the correlation observed in simulation of supercooled liq- 
uids In particular, slow and fast dynamical hetero- 
geneities correlate with regions of low and high energy, 
respectively. Since each column of the pictures in Fig. 2 
is an equilibrium configuration of the noninter acting H, 
the structure seen in the trajectories is purely dynamical. 

Down spins must form closed domains in a trajectory 
as a consequence of the local and causal nature of the dy- 
namical constraints. This fact is apparent from the illus- 
tration in Fig. 3. A spin is able to flip only if an appropri- 
ate neighboring spin is in the up state at the same time. 
It therefore will have flipped up previously and/or it will 
flip down later. As such, a well defined closed boundary 
must exist between regions of up and down spins. These 
boundaries are formed out of segments with shapes like 
those depicted on the top of Fig. 3, with the possibility 
of different slopes. In the FA model, the dynamical con- 
straint is spatially symmetric, and four possible kinds of 
boundary segments are allowed, with the restriction that 
all the segments of the first (second) kind must be be- 
low those of the third (fourth) kind. It follows that spin 
down domains must form semi convex polygons like that 
pictured at the bottom left of Fig. 3. Namely, any spatial 
line cuts the boundary only twice. In the case of the East 
model, only the third and fourth boundary segments are 
allowed, which restricts the domains to shapes like that 
pictured at the bottom right of Fig. 3. Since slow dy- 
namical heterogeneities correspond to spatial projections 
of the compact space-time down spin domains, they are 
necessarily compact, while the converse is true for fast 
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regions. This observation is in agreement with what is 
found with computer simulations of atomistic models j| . 

The geometrical construction described above implies 
that spin down domains cannot penetrate each other. 
Therefore, the set of trajectories maps to the configura- 
tion space of a two dimensional mixture of polydisperse 
non-interpenetrating objects of all the possible shapes 
allowed by the dynamical constraints. This observation 
motivates a description of the dynamics in terms of p(7, t), 
the density of domains of typical spatial size (height) / 
and extension in time (length) t. The partition function 
for trajectories is then Z ~ f Dp(l, t) exp (— Q[p(l, £)]), 
where Q[p(l, t)] is the free energy functional of this den- 
sity. Whatever estimate is used for ft [p (Z, £)], it is crucial 
that it enforces the constraint 

J dtp(l,t)=p eq (l)=ce- cl . (2) 

This condition ensures that any spatial cut of the equi- 
librium trajectories is an equilibrium configuration. For 
a density functional theory, it provides an approximation 
to the effects of the antiferromagnetic interactions of AS. 

In general, p(l,t) = p(l\t)p(t), where p(l\t) is the prob- 
ability density of I conditioned on t, and p(t) is the prob- 
ability density of domains of length t. In the case of the 
FA model, a domain is bounded by the random walks of 
two up spins between successive encounters, and p(l\t) 
can be obtained by standard arguments [fl5| 

p(l\t) ~ l 2 (Dt)~ 3/2 &cp(-l 2 /Dt) (3) 

with D ~ c. If the domains were isolated, p(t) would 
be given by the probability of first return of a random 
walker, which goes as p(t) ~ t~ 3 / 2 for large t, leading 
to (t) — > oo, and to the formation of unbounded do- 
mains. This result, however, corresponds to only the fer- 
romagnetic part of the interactions in trajectory space. 
It is frustrated by the condition ([|). By combining (^) 

with (§), we have instead p(t) - (rt) _1 / 2 e~V^, with 
D~ x c~ 2 - e 3 / T . This result for r is precisely the 
relaxation time in the FA model. The forms of p(l\t) and 
p(t) are verified in simulations of the FA model H]. We 
also see that in terms of the scaling variables I* = cl 
and t* = t/r, the density p(l*,t*) is independent of 
temperature. An illustration of these scaling relations 
is given in Fig. 4. The typical height of domains grows 
with the cubic root of the relaxation time, thus explain- 
ing the slow increase of the correlation length with 
decreasing temperature in the structure factors of Fig. 
1. Slow growth and domains being mesoscopic rather 
than macroscopic are consistent with experimental ob- 
servations J^]. Moreover, since the trajectories are ex- 
tensive in interfaces, the form of the structure factor for 
large k corresponds to a spatial projection of Porod's law 

In the case of the East model, typical distances 
and times are related at low temperatures by a T de- 
pendent dynamic exponent z(T) ~ l/(Tln2) |]l7| , |l2| . 



We therefore expect the conditional probabilities to be 
p(l\t) oc exp [—l z ( T ^/t]. The condition (||) then leads to a 
stretched exponential form for the persistence function, 
J t °° dt'p(t') ~ exp[-(£/r)^ T )], where (3{T) ~ l/z(T) 
and r ~ c~ z(yT \ in accordance with previous results [ [l8| . 
In this case, £a£ increases even more slowly with decreas- 
ing T than in the FA model. The geometry of the do- 
mains also explains the differences in the FA and East 
structure factors shown in Fig. 1. In contrast to the FA 
case, where neighboring domains can be compressed to 
the point of contact, the boundary between domains in 
the East model cannot be formed just by a single line of 
up spins in space-time. Instead, to decrease energy and 
thus increase probability, the boundary is wet by smaller 
domains, ideally like in Fig. 5. This structure has fractal 
dimension df = In 3/ In 2, and gives rise to the behavior 
S(k) ~ k~ df at intermediate momenta. At larger /c, S(k) 
probes the granular structure of the boundary, and thus 
exhibits the oscillatory behavior shown in Fig.l Jl9|| . 

The ideas presented here in detail for the FA and East 
models generalize to higher dimensions provided the dy- 
namical constraints are causal and local, so that excita- 
tions favor the creation of neighboring excitations. In 
general, dynamical constraints will appear only below a 
crossover temperature T x or above a corresponding pack- 
ing fraction. This can be incorporated into our descrip- 
tion by allowing the coupling A to be a function of tem- 
perature or packing fraction, increasing with decreasing 
T or increasing packing. The temperature T x will then 
coincide with the so-called "landscape" temperature [Q. 

The domains that thus appear for T < T x and grow 
with further decrease in T can lead to a broken sym- 
metry. Consider, for example, the partition function 
Zt^T) for trajectories between configurational fluctu- 
ations with overlap q at time difference T, Z^(<7, T) = 
E{a ,ar} P(°o) Z *o,Jt & {q ~ Ei ^o&T iT ), where 
San = (Jit — (cr), and p is Boltzmann distribution. 
Since Zt (<?, T) is proportional to the number of such 
trajectories, it is natural to write it as Zt^T) oc 
exp[N lut(q^ where ^(g, T) is the entropy density 
in trajectory space. It is approximately the entropy of 
mixing of the slow domains in trajectory space. At values 
of T for which the temporal extension of these domains 
is much smaller than T, lot (q, T) is extensive in time, 
and Zriq, T) is peaked at q = 0. In this case, correlation 
functions are exponential. 

Non-zero overlap is probable, i.e., Zt (q,T) is peaked 
at finite g, only when T is comparable to or smaller than 
the length of typical slow domains. As temperature is de- 
creased, the size of these domains increases, constricting 
the available trajectory space. At a low enough tem- 
perature, we may therefore expect out (q, T) to become 
sub-extensive, leading to a probable finite q throughout 
a relatively large range of T. This change in behavior 
of ujt (<?, T) signals a dramatic change in the number of 
available trajectories and a corresponding onset of dy- 
namical arrest. In this picture, therefore, the glass tran- 
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sit ion coincides with an entropy crisis in trajectory space, 
rather than in configuration space. 
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The idea of wet interfaces in configuration space appears 
in a recent equilibrium theory of glass formers Q , but it 
is difficult to see a connection between that idea pertain- 
ing to statics and the wet interfaces of trajectory space 
discussed here. 
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1. Structure factor S(k) of sf(t;At) for the FA 
(left) and East model (right). The S(k) are the spatial 
Fourier transforms of the normalized correlation functions 
(s 2 i (t;At)s 2 j (t;At))/(st(t;At)) (j = i, i ± 1, . . .). (•) indicates 
equilibrium ensemble average, so S (k) is independent of t. 
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FIG. 2. Equilibrium trajectories at T — 1.0 for the un- 
constrained case (top), the FA (bottom left) and East model 
(bottom right). The vertical direction is space, corresponding 
to a spatial window of systems of size L = 10 5 . The horizontal 
direction is time. Black/ white correspond to up/down spins. 




FIG. 3. Geometry of slow domains imposed by the dynam- 
ical constraints. Top: allowed boundaries between regions of 
up (black) and down spins (white). Bottom: shape of domains 
in the FA model (left) and in the East model (right). 
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FIG. 4. Equilibrium trajectories in the FA model for 
T = 0.6, 0.5, 0.4. Vertical direction is space, which scales 
with 1/c, and horizontal is time (t/10 3 ), which scales with r. 




FIG. 5. Energetically favored wetting of the boundary of a 
spin down domain in the East model. 
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